A non-linear model of trading 
mechanism on a financial market 

N.Vvedenskaya 1 , Y.Suhov 1,2 ' 3 , V.Belitsky 3 

Abstract 

We introduce a prototype model in an attempt to capture some as- 
pects of market dynamics simulating a trading mechanism. The model 
description starts with a discrete-space, contiinuous-time Markov process 
describing arrival and movement of orders with different prices. We then 
perform a re-scaling procedure leading to a deterministic dynamical sys- 
tem controlled by non-linear odinary differential equations (ODEs). This 
allows us to introduce approximations for the equilibrium distribution of 
the model represented by fixed points of deterministic dynamics. 

1 Introduction 

This paper proposes a model that takes into account, in a rather stylized form, 
some aspects of automated trading mechanisms adopted in modern financial 
markets, in particular, the dynamics of the limit order book. 

In short, a limit order book keeps records of arivals, movements and depar- 
tures of market participants (traders) who declare their trading positions. An 
arriving trader may wish to buy or sell at a certain price, and can move his 
declared price when time progresses. If the declared price is met by a trader 
with the opposite intention, a trade is recorded: this may lead to disappearance 
of one or both participants from the market, due to exhaustion of their offers. 
For a detailed description of some common limit order book models and their 
applications, see [3j 0] and references therein. 
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In the current paper we present a somewhat different model, including el- 
ements of queueing behavior of arrived offers; this model is studied by using 
techniques of asymptotic analysis. An earlier account of this work (in its pre- 
liminary version) can be found in [T]. Compared with 1 , in the present text 
we adopt a continuous-time setting for the basic Markov process: it clarifies the 
meaning of the main parameters of the model and shortens the proof of some 
of our main results. 

The model under consideration is a prototype; at this stage it does not aim 
to take into account all possible aspects that can be viewed as defining, either 
theoretically or practically. Instead, we opted for a simplified description which 
leads to some straightforward, yet instructive, answers. 

Our model differs from known models of the limit order book in a number of 
aspects. Arguably, it can be a subject of criticism (as a number of other proposed 
models). In particular, the strategic behavior of the model in its current form 
(and some further details) do not quite match existing mechanisms governing 
electronic trading on financial markets. Nevertheless, the model shows a certain 
amount of flexibility, and covers a borad range of situations. Its mathematical 
advantage is that in the scaling limit under consideration, it leads to a single 
fixed point. 

Our scaling limit is based on suppositions that (i) the number of market 
participants is large, (ii) during a very short time period only part of them 
makes a decision of performing a trade or maiking a move along the price range, 
and (iii) the probability for any given participant to make such decision is small. 
This makes it natural to change, in a suitable manner, parameters of original 
Markov process. 

After rescaling, a limiting dynamical system emerges, with a determinis- 
tic behavior described by a system of non-linear ordinary differential equa- 
tions. The rescaling techniques greatly simplify the structure of the model, 
and this phenomenon extends far beyond basic examples like the current proto- 
type model. As we mentioned earlier, the present paper focuses on a simplified 
model, with 'minimal' number of constant parameters, where some of technically 
involved issues are absent. 

A similar approach is commonly used in the literature on stochastic com- 
munication networks; see, e.g., [5j [6] [7] [8] and [9]. We also find similarities, 
as well as differences, with models proposed (in a different context) in a recent 
paper [10]; analogies with [10] could be useful for the aforementioned purpose 
of defining the prices that are appropriate for trades. 

In the next section we describe the underlying Markov process. In Section 
3 the rescaling of the process is presented and the main results are stated and 
the proofs are given. The last section contains concluding discussion of various 
aspects of the model. 
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2 The underlying Markov process 



The rationale for the models below is as follows. We consider a single-commodity 
market where prices may be at one of N distinct levels (say, c\ < ci < . . . < cjv, 
although the exact meaning of these values is of no importance here) . 

The market is operating in continuous time t £ K + where R+ = [0, oo). (As 
was mentioned above, the earlier version [T] used a more cumbersome discrtete- 
time version of the underlying process.) At a given time t £ M+, there are bi{t) 
traders prepared to buy a unit of the commodity at price Ci and Si(t) traders 
prepared to sell it at this price, which leads to vectors 

b{t) = (bx(t),...,b N (t)), s(t) = (8i(t),...,8 N (t)) £ Z%. 

Here and below Z + = {0, 1, . . .} stands for a non- negative integer half-lattice 
and Z+ for the non-negative integer A-dimensional lattice orthant. The pair 
(b(t), s(t)) represents a state of a Markov process {U(t)} that will be the subject 
of our analysis. 

Suppose that, for given k = 1,...N and t £ K + , we have that bk(t) > 
Sfe(t) > then each of the sellers gets a trade at a given rate px > and leaves 
the market, together with one of the buyers. Therefore, both values bkit) and 
Sfc(t) decrease at rate px- In addition, each one among Sk{t) sellers (i) quits 
the market at rate pq > or (ii) moves to the price level Ck-i at rate pm > 0, 
if k > 1. Similarly, every buyer among the bk{t) buyers (i) quits the market at 
the same rate pq as above or (ii) moves to the price level Ck+i, again at rate 
Pm, provided that k < N. 

Symmetrically, if Sfc(i) > bf.(t) > then each of the buyers gets a trade 
at rate px and leaves the market, together with his seller companion. The 
remaining traders at the price level proceed in a manner as above. 

Further, when fc = N , a buyer leaves the system with rate pq + pm- Simi- 
larly, for k = 1, a seller leaves the system with rate pq + pm- 

Finally, a random Poisson flow of new exogenous buyers arrives at the price 
level ci; the rate of this arrival equals Ab > 0. Similarly, a Poisson random flow 
of new sellers arrives at the price level cat; the rate of this arrival is As > 0. 

As usually, standard independence assumptions are in place. This generates 
the aforementioned Markov process {£/(£)} with trajectories {(b(t), s(t))\, t £ 
K+. 

Theorem 1 For any values of parameters Ab/s? Pq/m an d Pt, the process 
{U(t)} is irreducible, aperiodic and positive recurrent. Therefore, it has a unique 

set of equilibrium probabilities tt = (jr(b, s) : b,s £ Z^^j, and for any initial 

state U(0) (deterministic or random), the distribution of the random state U(t) 
at time t converges weakly to tt as t oo: 

lim P(z7(t) = (b,s)) = 7r(6,s). 
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Proof of Theorem^ Irreducibility and aperiodicity of the process is evident. 
Positive recurrence follows from the following observation. The (random) time 
that a given trader (a buyer or a seller) spends in the system, i.e., the time from 
his arrival till exit, is majorized by a sum of N independent exponential vari- 
ables. Therefore, the process {U(t)} can be majorized, in a natural fashion, by 
an M/M/oo queueing process. But the latter is known to be positive recurrent. 

The remaining assertions of Theorem Q] are standard. A 

Despite a concise description, the detailed pattern of behavior of process 
{U(t)} is rather complex, particularly for large values of N. For instance, 
consider the differences between vectors bit!) and b(t) and between s(t') and 
s(t), on a time interval (t,t') where < t < t'. The increments for the entries 
6fe( • ) and Sk( • ) for 1 < k < N — 1 are captured by the following equations: 

bk(t') = b k (t) + i^it, t') - n k (t, t') - i fc M (i, - ifc Q (;, 0, 

and 

s k {t') = *fc(*)+j^ 1 (t,0-»fc(*,0-j* M (*,*')-i fc Q (*,*'). 

Here ijf'it, t') is the number of buyers who move within time interval (t, t') from 
level k to k + 1 and j^it,t') that of sellers who move from level k to k — 1. 
Next, 4') is the number of buyers who quit the system during interval (t, t') 
from level and j^{t,t') the number of sellers who quit the system from level 
k. Finally, rik(t,t') is the number of buyers and sellers who got a trade over 
(t, t') at level k. All listed quantities are non-negative integer-values random 
variables. For k = 1 the structure of the expression is similar, with the term 
i^i(M') being replaced by i B (*,*') > 0, while for k = N the term j fc + x (*,*') is 
replaced by js(M') > 0; both ib(M') an d js{t,t') being distributed according 
to a Poisson law with mean Ab/s(^' ~ t)- 

In the simplest case of a market with one price level (N = 1), the pro- 
cess {U(t)} is a continuous-time random walk on the two-dimensional lattice 
quadrant Z^_, where 

bit') = bit) + i B it,t')-nit,t')-i Q it,t') 

and 

sit') = sit) + j s (M') ~ nit,t') - j Q it,t'). 

This already makes analytical representations for the invariant distribution n 
rather complicated; cf. QT] and references therein. 

3 Scaling limit 

The complexity of the time-dynamics and of the equilibrium distribution n for 
process {[/(£)} makes it desirable to develop efficient methods of approximation. 
In this paper we focus on one such method based on scaling the parameters of 
the process (including states and time-steps). 
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The re-scaling procedure is as follows: we fix values 7 > 0, /? > 0, a > 0, 
Ab > and As > and set: 

7 j3 ct 

PT = p Q = -, p M = y- (1) 

In addition, we re-scale the states and the time: pictorially, 

h s k t 

x k ~j-, y k ~j-, t~-. (2) 

Formally, denoting the Markov process generated for a given L by U^ L \ we 
consider the continuous-time process 

V (l \t) = \u {l) (tL), t>0. (3) 

Let denote a positive orthant in N dimensions. Suppose we are given 
a pair of vectors (a;(0),y(0)) e x where x(0) = (xi(0), . . . , x N (Q)), 
y(0) = (2/1(0), J/jv(0)). Consider the following system of first-order ODEs 
for functions x k = x k {r) and y k = 2/fc(r) where r > and 1 < i < N : 



±1 = A B - [0 + a) xi -7 min [a;i,j/ij, 
i fe = aifc-i - \B + ajx k - 7 mm [x k ,y k ], 1 < fc < N, 



(4) 



2/fe = ayfe+i - \ P + aj Vk - 7 min [xk, Vk] , I <k < N, 
y N = As - (p + o^ypf -7min[^jv,2/Ar], 

with the initial data 

x fc (0) > 0, y k (0) > 0, 1 < k < N. (5) 
The fixed point (x*,y*) of system has 

a;* = (a£,...,a;^) and y* = (y x , . . . , y* N ) 
where x* k and y k give a solution to 

Ab = \d + a\xl +7 min [a:*,y*], 

1 = [p + a) x* k + 7 min [x* k , yj] , 1< fc < JV, < 2 ' > 

"2/fc+i = + + 7 min [x* k ,y* k ], 1 < k < N, ^ 

As = [P + a)y* N +7min [x* N ,ylf\. ( - 4 "> 



(6) 



(In © we noted individual equations by addition signs that are used below). 
Both systems Q and © are non-linear. However, the non-linearity 'disappears' 
at a local level which greatly simplifies the analysis of these systems. 

In Theorems 2 and 3 below, we use the distance generated by the Euclidean 
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Theorem 2 (a) For any initial data (x(0),y(0)) with 2^(0) > 0, j/i(0) > 0, 
1 < i < N for all r > there exists a unique solution (cc(r), y(r)), to problem 
©-© and Xi(r) > 0,y;(r) > 0. 

(b) j4,s t — > oo, the solution approaches a fixed point which is a unique 
solution to system 



dist 



x(r),y(T)),(x*,y*) -+ 0. (7) 



Proof of Theorem® (a). Obviously, a unique solution exists for sufficiently 
small r > 0. Suppose that = max Xk and max Xk < Xi. Then xi < 0. 

k k<i 

Similarly if yj = max yk and max yk < yj then yj < 0. Also, Xi < if 

£i > Ab/(/3 + a) and t/at < if y^ > As/(/3 + a). Moreover, ±\ > when 
£i = and yN > when yjv = 0. Thus the components of the solution Xiir) 
and yi(r) are non- negative and uniformly bounded. By standard constructions 

of the ODE theory, a unique solution |(a;(r), y(r))| exists for all r > 0, and 

z(r),»(r)6B^. A 

Before proving assertion (b) , we discuss several properties of the solution to 
(SJ) - (0) . The following Proposition [1] indicates that the solution to ((4]) possesses 
a kind of the min/max principle. 



Proposition 1 Suppose that we have two initial points, x — (xi, . . . ,x n ) and 
x' = (x[, . . . , x' n ) £ such that for some k = 1, . . . , N, 

Xk < x' k , and y k > y' k . (8) 

Then for the solutions (x(t), y(r)), (x'(t), y'(r)), to (j4]) 7 with the initial con- 
ditions (as(0),y(0)) = (x,y) and (x'(0), y'(0)) = f or al1 T > °> 

Xk(r) < x' k (r), and y k {r) > y' k (r). (9) 

Similarly if Xk > x' k and yk < y' k then the corresponding solutions obey Xk (t) > 
x' k (r) and y fc (r) < y^if), for all r > 

Proof of Proposition [TJ It is sufficient to consider the situation where the 
strict inequalities take place. Suppose that (J9)) holds strictly for t < tq and 
fails at t = tq. For instance, let Xk(ro) = x k (To) and assume that k > 1 is the 
minimal index for which such an equality takes place. 

If yfe( T o) > y'k( T o) then Xfe(ro) < x' k (ro) and the above equality Xk(ro) = 
x'Atq) is impossible. The other case is considered in a similar manner. ▲ 

A corollary of proposition [T] is 

Proposition 2 If x k (0) > and yk(0) < for all 1 < k < N then x k {r) 
increases and J/fc(r) decreases in t, for all r > and 1 < fc < N. Similarly if 
ifc(O) < and yk(0) > 0, 1 < k < N, then Xk(r) decreases and J/fc(r) increases 
in t. 
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Proof of Proposition® It again suffices to assume that the strict inequalities 
hold true: Xfc(O) > and y/c(0) < 0. Then, for a small S > 0: x k {8) > x k (0) and 
Vk{8) < y k {0). 

Set x' k (0) — Xk(5), y' k (0) = Vk{8)- The coefficients of equations do not 
depend on r, therefore x' k {r) = x k {r + S) and y' k {r) = j/fe(r + S) for all r > 0. 
By Proposition 1, x k (r + S) = x k (r) > x k (r) and y' k (r + S) = y' k (r) < y k {r), 
for all r > and 1 < k < N. As 5 may be arbitrarily small, the assertion of 
Proposition [5] is valid. ▲ 

Proof of Theorem 2, (b). Given x(0) = (xi(0), . . . ,x N (0)) and 
2/(0) = (yi(0),...,2/iv(0)) € K+ , let (x(r),a5(r)) be the solution to ©, ©. 
Consider two additional solutions, (cc'(t), y'(r)) and (x" (t) , y" (t)) , to Q with 

= °, yfc(°) = max A s /(a + maxyi 
and 

4'(0) = max \\ B /(a + p), max a* (0)1 , y k (0) = 0. 

L i J 

By Propositions [1] and [5] 

zUr) < x fc (r) < 4'(r), ^(r) > j/*(r) > y£(r)- (10) 

Further, xi(r) increases, while y' k (j) decreases in t. By the same token, x' k (r) 
decreases and y k (r) increases in r. Therefore, both pairs (x'(t), y'(r)) and 
(x"(t), y"(r)) tend to limits as r — s> oo, which are fixed points, i.e., solutions to 
Eqns ©■ By (ITU1) . any solution to ((3]), ([S]) eventually lies between these limits. 
To finish the proof of the theorem, we have to show that the solution (x* , y*) 
to Eqn © is unique. 

For convenience, we state the corresponding assertion as Lemma [1] 

Lemma 1 For any values Ab/SjQ^T and ft > there exists a unique solution 
to Eqns (0|). 

Proof of Lemma\\$ To start with, note that every solution to © has 

i* > x* 2 > . . . > x* N , yl < 2/2 < . . . < y* N . 

Therefore, if x\ < y\ then x k < y k , 1 < k > N and if xn > y^ then x k > 
y k , 1<> k < N . 

It is convenient to introduce auxiliary variables v k ,w k > in terms of which 
Eqns ([5]) will be treated. Geometrically the idea is as follows: we start with 
Eqn @W : A B = {a + P)v x + min [vi, wi] and watch how this relation between 
v\, w\ is transformed by ([B])^^,®^ 3 ^ to relation between i> 2 , u>2, then to relation 
between U3, W3 etc. 

The locus of points (it, v) S R 2 where 

U, u> > 0, Ab = (a + j3)v + min[u, w] 
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coincides with a continuous broken line, L\ C M. 2 ,, formed by two pieces: 1) a 
vertical ray Lf^ emitted from the point lying on the bisectrix where 

o[ 0) = (X B /(a + f3 + 7), A B /(a + /3 + 7)) 



and 2) a line segment of a negative slope — (a + /3)/7, joining with the 
point on the horizontal axis: 



(o) 



«1 



(A B /(a + /3),0). 



See the figure below. 




V«*+P) 



The passage from values v\, w\ to «2, W2 generates a 1 — 1 map acting on 
L\. The image of L\ under the map {y\,w\) i-> («2,W2) is another continuous 
broken line, L 2 C M+, formed by three pieces: 1) a vertical ray issued from 
point 



(X B a/ (a + 13 + 7) 2 , A B /a), 



2) a line segment joining the points and ckp where lies on the 
bisectrix and has both co-ordinates equal to 

X B a(a + fi + 7) 



(a + (3 + 7) [(a + f3)(a + j3 + 7 ) 2 + a 2 7 ] 



and 3) a line segment joining o^ 1 ' with the point a^' on the horizontal axis: 



(2) 



J 2 ) 



\ B a 



. 
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The slopes dw/dv of segments L% ' and l\ are negative, but the slope flattens 
when we pass from to l\ ■ 

In the above figure, L\ and L2 are shown on the same (u,w)-plane R 2 ; in 
this figure line Li lies to the left of L\. (The third broken line, Mm present in 
the figure is explained below.) 

A similar picture persists when we iterate, i.e., pass from (u 2 , Wi) to (i>3, w^) 
and so on. At step k {vk,w^) i-> (Vk+i, Wfe+i) where again the map is 1 - 1, 
(ufc+i, Wk+i) belongs to a continuous broken line Lj~ + i in R 2 ^ formed by k + 2 
pieces. One piece, i^+i' i s a vertical ray while the k + 1 others, L^_^, . . ., 
, are line segments of negative slopes, the slope flattens when we pass 

from iJjp to L^. The last segment, £&i , joins a point on the bisectrix and 
a point on the horizontal axis. 

At the end of this process we obtain a continuous broken line Ln, the locus 
of points (vn,wn)- Our next step is to consider the intersection of Lm and Mm 
where Mm is the locus where 



More precisely, Mm is a continuous broken line formed by a horizontal ray issued 
from the point 



lying on the vertical axis. Cf. the figure. 

We want to check that the point of intersection is always unique: it yields 
to unique solution to ©. 

Line Lm may intersect the horizontal part of Mm- That means that there 
exist a solution to ([6]) where x\ > y\ and x* N > V*m (this case is not presented 
on our figure) . For the proof of lemma it is needed to show that in this case Lm 
cannot intersect the sloppy part of Mm, where dm > xm ■ For that it is sufficient 
to prove that dw /dv on Lm \ L^ i t ne part of Lm above the bisectrix, is steeper 
than —7/ (a + /?), the slope of the segment of line Mm- 

On the other hand if Lm does not intersect the horizontal part of M m it has 
to intersect the sloppy part of Mm and it is needed to show that in this case 
such intersection is unique. Here again it is sufficient to show that dw/dv on 

Lm \ L^ is steeper than —7/(0 + ft). 

To prove the assertion of the lemma we show that dw/dv on Lm \ L^ is 
always steeper than —7/ (a + (3). In fact, it suffices to verify that 



As = {a + (3)wn + 7 mm [vm, wm]- 



Va + /? + 7' a + /3 + 7'/ 
and a line segment joining this point with the point 





dw m 
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on L N , l<i< N. 



(11) 



dv N 



< 



a + p 
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Any segment L k % \ 1 < i < k, maps onto segment segment L k k ' maps 

onto two segments: L k k \ and . The slope of segment L^\ x ,i = 1, . . . , k+1 

is steeper than that of L k 1 ' and the slope of L^+P is steeper than that of L k k K 

(i) 

More, the slopes of L y N , < i < N, are steeper then the steep segment of Mjv- 
To show that consider three cases of the map Lk — > Lk+t- 

1) Vk > w k , v k+ i > w k +i; 

2) v k > w k , v k +i < w k+1 ; 

3) v k < Wk+x, v k +i < Wk+i- 

In these cases we have, the following equations 

( A k+1 (a + /3 + ~f) dw k 



dwk+i 



dv 



fe+i 



- = I (a + 13 + if dw k 

dv k ' 



dv k 



-B fc (q+/3+7) dw k 
a 2 dvk 



(12) 



where, respectively, 



(a + /3)v k+ i +jw k+ i (a + (3 + i)w k 
Vk = , Wk+l = : 

(a + (3 + i)v k+1 (a + (3 + j)w k 

Vk = , w k+ i - 



Vk 



(a + f3 + j)v k+1 



a 



Wk+l 



{a + P)i 



Wk + iv k 



Here A k+ i = a + (3 + 7 k+1 and B k — a + (3 + 7—. 

v k+ \ w k 

Therefore for N > 2 



dw 



N 



dv 



N 



{a + 13 + 7) JV " 1 dtoi _(a + /3 + 7) JV " 1 a + 



vN-l 



dw 



dv\ 



JV-1 



dw 
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a + (3 



For iV = 2 we have the middle case in (1121) with k 
lin the i 
we get that for N = 2, 



again the needed inequality takes place. In fact, 



: 1, k 

7 



+ 1 = 2. Here 
By using (|T3 



(a + /? + 7) 2 dwi 



on L3,. 



This finishes the proof of Lemma [TJ A 



Theorem 3 Suppose that the re-scaled initial states converge in probability: for 
any e > 0, 

It/(0M,r(0)..y(())) > < | = 0. 



lim P dist 
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Then, for all T > 0, the process ^V^ L \t), t £ [0,T]| converges in probability 
to the solution |(x(r),j/(r)), < r < t|. T/iat is, V e > 0, 



lim P I sup dist 

L->oo \o< T <T ^ 



VW(t),(x(t),»(t))]} 
In particular, if x(0) = x* and y(0) = y* i/ien 



> e = 0. 



lim 



( sup {dist\v {L) (r),(x*,y*)]}>e) 

\0<T<T 11 J ' J 



> e = 0. 



(13) 



(14) 



L^oo \o<r<T 

Moreover, if process \ U(t), } is in equilibrium then Eqn |i^| ) ZioZds trite. 

Proof of Theorem [3j Let G^ L ' denote the generator of the Markov pro- 
cess {U^ L \t)} (with rates as in Eqn (pj). Then the action of matrix G^ L ' on 
functions 0(6, s) of state variables 6, s £ lJ^_ is determined by the equation 

G* £ V(M) = | E ([^(6-e fc ,«)-^(b,«)]6 fc 

l<fc<JV 

+ [<p(b,s-e k )-<p(b,s)] Sk j 
^2 [0( b - e fc + e fc+ i,s) - 0(6, s)]6 fe 



l<k<N 



+ [4>{b, s - e k +e fe _i) - 0(6, s)]s k 

Kk<N 

A B [0(6 + ex, a) - 0(6, a)] + A s [0(6, a + e N ) - 0(6, a)] 

2 



(15) 



+ 7 H [0(6-e fc ,a-e fe )-0(&,a)](&fcAs fe ). 



l<fc<7V 

Here e k , 1 < k < N, stands for the vector in whose components are all 0's 
except for the fcth one, equal to 1. 

In particular, we can take a function = 0' L ) of the form 0(6, a) = 
6 s \ 

where ip is a smooth function on R+. This choice agrees with the 



'•P 



V L 



spatial scaling x ~ 6/L, y ~ s/L m Eqn (0). Then Vs = (xi, . . . , Zjv) € 
and y = (y 1; . . . , y^r) G R+ j with 6 = [-^J , s = L-^yJ where [ • J stands for the 
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integer part, we obtain that 

G {L ^{[Lx\/L,[Ly\/L) 

13 \[p((b-e k )/L,s/L) - V {b/ L,s/ L^b k / L 

+ [ V (b/L, (s - e k )/L) - <p(b/L, s/L)] s k / 



Kk<N 



{ [p((b-e k + e k+1 )/L,s/l?) -<p(b/L,8/L)]b k /L 

^ l<k<N 



Y [<p(b/L,(s-e k + e k ^)/L) - ip(b/L,s/l?J s k /L 

Kk<N 



+A B 
+A S 



+ ei)/i, s/i) - s/i) 
(s + ejv ) /l) - <p(b/L, s/L 



+7 



5^ [^((b-e fc )/L,(s-e fe )/L) - <p(b/L, (6* A s k )/L . 



Kk<N 



Next, we multiply the both side by L - in agreement with the time-scale r ~ t/L 
in Eqn @ - and pass to the limit L — >• oo. This yields 

lim LG (L) (p{x,y) 

L—>oo 



[ -13 Y (x k —+y k — 

r / a 3\ 



+CK| ^ ^ 
"l<fc<AT 

+ f Abt^- + A s 



K^., Xk ^dx k+1 dx k ) ,?rt„\dy k -i dy k 



d 



■dx k+ 

a 



Kk<N 



dx\ dy 



l<fc<JV ° Xl ° Vl J 



Applying Theorem 6.1 from [T2], we obtain Eqn (fTU|). 

The next remark is that each scaled process has a unique invariant 

distribution 7?^; the family of probability distributions w^ L ' (considered on 
Rtj is compact in the sense of convergence in probability. This can be deduced 
from the above remark that the original processes {U^}, and hence, the scaled 
process can be majorized by suitable analogs of M/M/oo queueing sys- 

tems. It is easy to see that every limiting point for tp l ) when L — > oo is a 
delta-measure sitting at a fixed point for system (|5|). However, the latter is 
unique and coincides with (x*,y*). Consequently, the distributions con- 
verge in probability to the aforementioned delta-measure. Then, applying the 
already established result, we obtain Eqn (fl^jl . ▲ 
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4 Fixed points in the scaling limit. Concluding 
remarks 



The approximation developed in Theorem 3 calls for an analysis of solutions to 
©• 

The parameter space formed by 7, chq/mj an d Ab/ S is partitioned into 
open domains where one of the following generic patterns persists: 

(i) vn > wn, (h) vi < w%, and (iii) > Wi for i = 1,...,£ and Vi < Wi for 
i = £+ 1, . . . , N where 1 < £ < N. In each of these domains system (j6j is linear. 

A particular algorithm for calculating (a:*,y*) is based on the following 
recursion. Set 

x (0) , (0) 

(0) _ A fc (0) _ Q M^_i (0) _ Ag (0) _ OiMU, + 1 

"Q + «M + 7 a Q + a M + 7 + "M CtQ + a M 

Next, let (x( k \ y^ k '\ k = 1,2, . . . be the solution to the system 

\ ( , A (fe) , • r (fe- 1 ) ( fc -!)i 

(fe) / , \ (fe) , • r (fe- 1 ) (fe- 1 )! 1 ^ • ^ i\t 
a m xlJi = [uq + a m \x\ + 7mm [x\ ' ,y\ '\,l<t<N, 

OmVi+i = \Uq + amjl/f + 7min [x[ k \y < i k " 1) ], 1 < i < N, 

A s = (a q + a m j V { n + 7 min > 2/yv ■ 

These iterations converge because the inequalities a;^ > x l - fe ~ 1 ' 1 , y^' < 
y^ ^ hold true V i, k > 1 and, values x^\y^ are uniformly bounded and 
there exist lim , lim y^ that, naturally, satisfy ©• 

fe— 5-00 fe— fOO 

We conclude with the following remarks. 

Our model presents also a caricature of "overproduction crisis": In fact, if 
A s is sufficiently large, so that Xi < y i7 1 < i < N, then the amount of "trades" 
J2i 7 rnin[xi, yi] — J2i~f x i ^ s no ^ changing by increase of A s , all "extra" sellers 
leave the market without performing any trade. 

It is interesting to investigate the dependance of trade performance on 7 
as 7 — > 00, that is where the trade action happens almost immediately after 
the moment when traders appear at some price level. Then almost all trades 
happen at two levels io and io + 1, and Xi is very small as i > ig + 1, yi ia 
very small as i < io. But our limit model does not permit to consider the case 
7 = 00, though, sure the initial Markov process can be investigated in case of 
immediate trade deals. The limiting case 7 = 00 of our model is close to the 
problems investigated in [ID] . 

We hope that the variation of these models parameters can help to deter- 
mine factors attracting or repelling various 'market participants'. An important 
aspect of any model of the market is what possibilities it gives for an accurate 
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prediction of the stochastic component in the dynamics of the market prices 
and volumes. 

The current set-up of the model presented here admits straightforward gen- 
eralizations to the case where parameters 7 and czq/m depend on i, < i < N 
and on the trader type (b/s). Another generalization emerges if these parame- 
ters and Xb/s become state-dependent. It is also possible to allow the exogenous 
buyers and sellers to enter the system at any price level among c±, . . ., cn- To 
take into account elements of the FCFS discipline, one could introduce various 
priorities into the dynamics of process U{t). 

Finally, we would like to note that there are several forms of convergence for 
which the assertion in Theorem 3 holds true. 
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